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Abstract 

The EDELWEISS experiment has improved its sensitivity for the direct 
search for WIMP dark matter. In the recoil energy range relevant for WIMP 
masses below 10 TeV/c 2 , no nuclear recoils were observed in the fiducial 
volume of a heat-and-ionization cryogenic Ge detector operated in the low- 
background environment of the Laboratoire Souterrain de Modane in the 
Frejus Tunnel, during an effective exposure of 7.4 kg-d. This result is com- 
bined with the previous EDELWEISS data to derive a limit on the cross- 
section for spin-independent interaction of WIMPs and nucleons as a function 
of WIMP mass, using standard nuclear physics and astrophysical assump- 
tions. This limit excludes at more than 99.8%CL a WIMP candidate with 
a mass of 44 GeV/c 2 and a cross-section of 5.4xl0 -6 pb, as reported by 
the DAMA collaboration. A first sample of supersymmetric models are also 
excluded at 90%CL. 



Introduction 

The experimental efforts in the search for Cold Dark Matter in the form 
of Weakly Interacting Massive Particles (WIMPs) are steadily increasing (see 
e.g. Ref. [0 for a review). In direct searches, the experimental signature of 
the WIMPs from the galactic halo is the observation of nuclear recoils induced 
by their scattering. Current experimental sensitivities for the interaction rate 
of WIMPs are of the order of 1 per kilogram detector material and per day 
for the various experiments at the forefront of this search 

The experiment DAMA has reported an annual modulation signal in 
Nal detectors. This represents a challenge to other detecting methods to 
reach equivalent sensitivities, a standard procedure to compare different ex- 
periments having been laid out in Ref. ||. Two experiments j|, ||, both using 
cryogenic heat-and-ionization germanium detectors, were able to exclude at 
more than 90%CL the central value deduced by DAMA from its annual mod- 
ulation signal for the WIMP mass and its nucleon scattering cross-section 
(M w = 52 GeV/c 2 and a n = 7.2 x 10~ 6 pb, respectively). The CDMS ex- 
periment was the first to report a limit excluding this value. However, 
the operation of the detectors in a shallow site, with only 16 meters of water 
equivalent (m.w.e.) protection from cosmic rays, leads to a sizable back- 
ground of nuclear recoils from neutron scattering that requires a delicate 
procedure of background subtraction J/J. The EDELWEISS experiment ||, 
located in a 4800 m.w.e. deep underground site, was also able to reject 
that value without requiring any background subtraction. However, the en- 
ergy resolution in that experiment restricted the sensitivity to nuclear recoils 
above 30 keV and the accumulated run time was not sufficient to extend 
the sensitivity to the central value obtained by DAMA |§ when their annual 
modulation signal is combined with their exclusion limit from pulse shape 
discrimination in Nal @ (M w = 44 GeV/c 2 and a n = 5.4 x 10~ 6 pb). Be- 
yond this, an important step in these searches would be to reach out to the 
Supersymmetric model calculations predicting the largest cr„ values. 

Following the results obtained by EDELWEISS using a heat-and-ionization 
cryogenic Ge detector ||, three new detectors were put in operation. The aim 
was to improve our understanding of the performance of such detectors, and 
to extend the sensitivity to lower cross-sections. This letter presents the im- 
proved cross-section limit achieved using a detector with improved charge 
collection and energy resolution. 

Experimental Setup 

The experiment is located in the Laboratoire Souterrain de Modane in 
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the Frejus Tunnel under the French-Italian Alps, under a 4800 m.w.e. rock 
overburden. The experimental setup is described in ||; only the relevant 
modifications are discussed here. 

Three new 320 g cryogenic Ge detectors Ui~0| , each 70 mm in diameter 
and 20 mm in height, are operated simultaneously. Each one is equipped 
with a segmented electrode defining two regions, a central part and a guard 
ring. To improve their reliability, all electrical contacts with the electrodes 
are ultrasound-bonded instead of glued as in Ref. ||. 

Two of the detectors, labeled GeA19 and GeAUO, are very similar to 
the one used in the year 2000 runs (Ref. labeled GeA16. The third 
detector, GGAl, differs by the presence of a 60 nm hydrogenated amorphous 
Ge layer deposited under the 70 nm Al electrodes and on all exposed surfaces. 
This modification was done to test whether an amorphous layer can improve 
charge collection properties, as suggested by [JEJ . Bias voltage values between 
2 and 4 V are used. 

The size of the NTD heat sensors and the thermalization of the detectors 
are improved for a better sensitivity, in light of the previous experience with 
GeA16. As a result, it was possible to operate the detector at a reduced 
temperature of 17 mK (regulated to within 10 /iK). 

The data acquisition system has been upgraded to a design with fully 
numerical data flow and trigger. The signals from the 3 heat and 6 ionization 
channels are continuously digitized at respective rates of 2 and 200 ksample/s 
and sent to the data acquisition PC via an optical link. The ionization data 
are then filtered on-line using an Infinite Impulse Response (IIR) high-pass 
elliptic filter of 4th order, in order to remove most of the microphonics noise, 
below a frequency of 1200 Hz. The trigger is defined by requiring a minimum 
threshold on the absolute value of any of the filtered ionization channels. 

Detector Calibration 

The heat and ionization responses to 7 rays were calibrated using 57 Co 
and 60 Co sources. In 2000 M, the performance of the detector GeA16 was 
partly limited by a poor baseline resolution on both ionization and heat 
channels. While resolutions at 122 keV of the new detectors remained close 
to those of GeA16, the baseline resolutions were somewhat improved. The 
ionization baseline resolutions are all below 1.5 keV FWHM, and are 1.3, 0.5 
and 0.4 keV for the heat channels in GGAl, GeA19 and GeAUO, respectively. 
The resulting improvement is illustrated in Fig. [I], showing the low-energy 
spectra recorded in the three detectors in the low-background physics run. 
Here, the energy corresponds to the average of the ionization and heat signals, 
weighted by the square of their respective resolutions. The 8.98 and 10.37 
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keV lines from the decay of the cosmic-ray induced long-lived isotopes 65 Zn 
13 1 are clearly resolved in GeA19 and GeAUO, with a resolution 
of 0.6 keV FWHM. The resolution in GGA1 is only 1.2 keV FWHM, but 
the two-peaked structure can again be observed. The degraded resolution of 
GGAl relative to the two other detectors is due to a reduced NTD sensor 
volume (1.6 vs 5.6 mm 3 ), and an increased sensitivity to microphonics of the 
center electrode ionization channel. 

The threshold level of the ionization trigger was measured using two dif- 
ferent techniques. The first one consists in extracting, as in Ref. ||, the 
threshold value corresponding to an efficiency of 50% from a fit to the low- 
energy part of the Compton plateau recorded with a 60 Co 7-ray source. The 
second technique was made possible by the simultaneous operation of the 
three detectors with a 252 Cf neutron source. Neutron scattering induces a 
large number of coincidence events where at least two detectors are hit. The 
efficiency curve as a function of ionization energy in one detector is given 
by the ratio of the ionisation energy distributions obtained with and with- 
out that detector appearing in the trigger pattern, the reference population 
being all events where at least one other detector took part in the trigger. 
Both Compton and coincidence techniques give consistent ionization thresh- 
old measurements within 0.2 keV. For GGAl, the values corresponding to an 
efficiency of 50% are 3.7±0.2 and 3.5±0.1 keV, respectively. This represents 
a significant improvement compared to the performance of GeA16, where the 
corresponding values varied between 5.7 and 11 keV during the run and re- 
stricted the analysis to nuclear recoils above 30 keV. With a 50% efficiency 
reached at 4 keV, the nuclear recoil selection described below reaches its full 
efficiency within less than 1% for recoil energies above 20 keV. 

The study of the distribution of the quenching factors Q (the ratio of the 
ionization signal to the recoil energy, calculated as in M) recorded in the 
presence of a 7-ray source revealed problems with the charge collection in 
GeA19 and GeAUO. Fig. ^| shows the Q distribution for recoils between 20 
and 200 keV recorded with 57,60 Co sources, for the different detectors. GeA19 
and GeAUO have been added together, since they display a very similar 
behavior. The distributions are normalized to the number of entries in GGAl. 
The distribution for GGAl is centered on 1, as expected by construction for 
electron recoils. The distribution for GeA19 and GeAUO show a narrower 
peak centered at 1, as expected from the better heat resolution. However, a 
relatively flat tail of events with Q values ranging from to 1 is observed. 
While 1.3 and 2.2% of events have Q values below 0.5 in GeA19 and GeAUO, 
this fraction is approximately 0.01% in GGAl. The tail amplitudes do not 
depend on the recoil energy range. On the basis of the 57,60 Co calibration 
runs it can be expected - and later experimentally verified - that these will 
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produce fake nuclear recoil events (Q ~ 0.3, see below) at a rate of a few 
events per kg-d in GeA19 and GeAUO. It was therefore decided that only the 
GGA1 data would be used for deriving a limit on WIMP interactions. 

Although the suppression of charge collection problems in GGA1 may 



indicate that the amorphous layer helps prevent them, as suggested by |12 
more thorough tests are needed before reaching any firm conclusion. 



Fiducial volume and acceptance 

As in Ref. [|S| , a fiducial volume is defined in order to exclude events occur- 
ring in the outer perimeter of the detector as it is more exposed to external 
sources of radioactivity and to charge collection problems. The selection cut 
is the same: more than 75% of the total charge must be collected on the 
center electrode. Here also, two methods are used to measure the fraction of 
the total detector volume thus defined. The first one uses the data collected 
with the 252 Cf source: the fraction of nuclear recoil candidates passing the 
fiducial cut is compared with the results of a Monte Carlo simulation of the 
neutron interactions in the detector. The second method exploits the unifor- 
mity of the 65 Zn and 68 Ge decays within the detector volume. The fraction 
of the total intensity of the 8.98 and 10.37 keV peaks selected by the fiducial 
acceptance is then equal to the fiducial volume fraction. The two methods 
give identical results for GGAl (57 ± 3%) and agree within 3% for the other 
two detectors. 

The acceptance for nuclear recoils is defined both in terms of ranges in Q 
and recoil energies. The neutron calibration of the three detectors confirms 
the parametrization used in Ref. ||, namely, the center of the band is given 
by O.IG(Er) 018 , where Er is the recoil energy in keV, and its width is equal 
to that predicted from the propagation of the heat and ionization resolutions 
added in quadrature with a constant rms spread of ~0.035. Again, the width 
of the band is set to ±1.645cr. It was verified on the neutron data that this 
selection does correspond to an efficiency of 90%. 

The lower bound of the recoil energy range for the selection of nuclear 
recoil is set to 20 keV, based on the same arguments as in Ref. ||: the ef- 
ficiency to nuclear recoils should be as uniform as possible within the band, 
and it should exclude regions where the 7-ray rejection, estimated by prop- 
agating the experimental heat and ionization resolutions, is expected to be 
worse than 99.9%. 

Given that a background of events with improper charge collection has 
appeared in GeA19 and GeAUO, and that it has a flat distribution in both 
Q and recoil energy, one could expect a similar behavior, albeit at a lower 
level, in GGAl. Therefore, some care must be taken in the definition of the 
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upper bound of the recoil energy range of WIMP candidates. The natural 
choice is to calculate, using the prescription of Ref. ||, the upper bound 
corresponding to 95% of all WIMP-induced recoils above 20 keV. Using the 
standard halo and nuclear form factor parameters^, this bound depends on 
the mass of the WIMP and varies from 33 keV at 20 GeV/c 2 to 86 keV at 100 
GeV/c 2 and saturates slightly above 110 keV at masses above 10 TeV/c 2 . 

Within the fiducial volume (57% of 318.5 g), the acceptance for nuclear 
recoils from WIMP interactions thus corresponds to 90% (width in Q) times 
95% of all recoils above 20 keV (mass-dependent recoil energy range). 



Results and Discussion 



The low-background physics data consists of all physics runs recorded 
over a period from February to May 2002. The physics data-taking period 
started a few months after the installation of the detectors in the Laboratoire 
Souterrain de Modane, after a period of optimization of resolution, calibra- 
tion and of long exposure to an intense 60 Co source while all electrodes were 
shorted. The running conditions were kept as homogeneous as possible, until 
the run was interrupted by an accidental warm-up of the detector. In addi- 
tion to a constant monitoring of the data, the homogeneity of the running 
conditions was checked with ^weekly 57 Co calibrations and two neutron cal- 
ibrations. The total physics run time at low background corresponds to 54 
days, of which 2.0% are lost due to the regular shorting of the electrodes to 
prevent the accumulation of space charge, 3.9% are lost due to the dead-time 
of the data acquisition and 6.2% are lost in a few hour-long episodes where 
the microphonics noise reached unacceptable levels, as attested by a strong 
deterioration of the baseline resolutions. The total exposure is thus 15.1 kg-d, 
of which 8.6 is in the center fiducial volume. The exposure corrected for the 
acceptance of the nuclear recoil band is 7.4 kg-d. 

The data recorded in the fiducial volume of GGAl are shown in Fig. |3|. 
The ionization-to-recoil energy ratios are plotted as a function of the recoil 
energy. Only events triggered by GGAl alone and with an ionization energy 
above 3.5 keV (hyperbolic dashed line) are shown. The 99.9% acceptance 
band for photons shown as a dotted line is the result of a simple propagation 
of the average heat and ionization resolution, assuming a Gaussian dispersion. 
The population of events around Q~0.5 associated with low-energy j3 and 



7 surface events, so prominent in our previous data of Ref. [TTfl and less so 



1 The halo parameters are a local WIMP density of 0.3 GeV/c 2 /cm 3 , a rms velocity of 
270 km/s, an escape velocity of 650 km/s and an Earth-halo relative velocity of 230 km/s. 
The Helm parametrization of the nuclear form factor is used with the recommended values 
of a = 0.52 fm, s = 0.9 fm and c = 1.23A 1 / 3 _ Q.6 fm. See Ref. M for details. 
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in Ref. ||, is only represented here by 4 to possibly 5 events. The correct 
interpretation of these events would require a significant increase in exposure, 
given the low level of background reached with the present detector. The 
same is true for the interpretation of the three events below Q=0.7 and with 
recoil energies between 119 and 182 keV. 

The event at 119 keV and Q=0.3 is lying at -1.646er of the centre of the 
nuclear recoil band. Given the uncertainty in the experimental determination 
of a, we conservatively choose to consider this event as a nuclear recoil, 
entering in the acceptance for WIMP masses above 10 TeV/c 2 . For lower 
WIMP masses, no events are in the nuclear recoil band. 

The absence of events in the defined acceptance is interpreted in terms of 
upper limits at 90% CL on WIMP-nucleon scattering cross sections for Mw 
< 10 TeV/c 2 following the prescriptions of Ref. || with the standard halo 
and nuclear models described above. 

The limit as a function of WIMP mass is shown in Fig. |], where it is 
compared with the previous result obtained with an effective exposure of 4.3 
kg-d of the GeA16 detector ||. The limit resulting from the combination of 
the two measurements, corresponding to an effective exposure of 11.7 kg-d is 
also shown. The 3a contour corresponding to the annual modulation effect 
of DAMA NaIl-4 || is shown: the black circle marks the central value of 
that measurement at M w = 52 GeV/c 2 and a n = 7.2 xl0~ 6 pb. The present 
combined results are incompatible with the interpretation of the modulation 
effect in terms of a WIMP behaving according to the standard phenomeno- 
logical model of Ref. 0. While 9.8 nuclear recoils should have been observed 
between 20 and 64 keV, none are observed. The Poisson probability of such 
a fluctuation is 0.006%. The black triangle on Fig. g at M w =44 GeV/c 2 and 
a n = 5.4 xlO -6 pb corresponds to the most likely value quoted by DAMA 
when they combine their modulation results with their limit achieved using 
pulse shape discrimination in Nal @. It is also incompatible with the present 
EDELWEISS results, the Poisson probability of observing no events from a 
prediction of 6.2 events being 0.2%. The only remaining part of the 3a NaIl-4 
DAMA zone corresponds to neutralino masses below the limit of 45 GeV/c 2 
obtained at LEP |14]| . Clearly, the standard prescriptions of Ref. Q fail at 
reconciling the EDELWEISS and DAMA experimental results. 

In Fig. |5| the combined EDELWEISS limit is compared to those obtained 
by other direct WIMP searches. The EDELWEISS sensitivity for spin- 
independent WIMP-nucleon interaction is the best achieved so far by any 
dark matter search for masses above 35 GeV/c 2 . Furthermore, the EDEL- 
WEISS data start to probe some of the supersymmetric models predicting 
the highest interaction rates. Fig. [| shows as an example the range of masses 
and cross-sections allowed in the calculations of Refs. ]T5L [Tj| , where relaxed 
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conditions of unification at the GUT scale yield higher upper bounds for o n 
(for comparison, see e.g. Ref. fL7|). 

Conclusion 

The EDELWEISS collaboration has searched for nuclear recoils due to the 
scattering of WIMP dark matter using a 320 g heat-and-ionization Ge detec- 
tor operated in a low-background environment in the Laboratoire Souterrain 
de Modane. After a combination with our previous data ||, the achieved 
sensitivity is so far the best for all direct searches for WIMP masses above 
35 GeV/c 2 . The limit obtained on WIMP-nucleon interaction cross-sections 
as a function of WIMP mass is based on the absence of events in the re- 
coil energy range relevant for WIMP masses below 10 TeV/c 2 and does not 
rely on any background subtraction. The combined EDELWEISS result is 
incompatible at more than 99.8% CL with a WIMP of mass 44 GeV/c 2 and 
a nucleon scattering cross-section of 5.4xl0 _6 pb reported by the experiment 
DAMA || based on the same standard nuclear physics and astrophysical 
assumptions. Furthermore, the EDELWEISS experiment excludes a first 
sample of supersymmetric models predicting the highest WIMP-nucleon in- 



teraction rates 15, 16 . 
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Figure 1: Energy pulse height spectra for low-energy gammas (sum of the 
ionization and heat channels, weighted by their resolution squared) in the 
fiducial volume of the EDELWEISS detector, for the low-background physics 
runs: (a) sum of the distributions of the detectors GeA19 and GeAUO; (b) 
distribution in the detector GGA1. The arrows indicate the peaks at 8.98 and 
10.37 keV, corresponding to the de-excitation of the cosmogenic activation of 
65 Zn and 68 Ge in the detector, and the 71 Ge activation that follows neutron 
calibrations. 
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Figure 2: Distribution of the ratio of the ionization pulse height to the recoil 
energy (quenching factor Q) obtained by exposing the detectors to 57 Co and 
60 Co 7-ray sources. The ratio is normalized to 1 for electron recoils using 
the photopeaks of the 57 Co source. Shaded histogram: detector GGAl. Line 
histogram: sum of the detector GeA19 and GeAllO, normalized to the same 
number of entries of GGAl. 
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Figure 3: Distribution of the quenching factor (ratio of the ionization signal 
to the recoil energy) as a function of the recoil energy from the data collected 
in the center fiducial volume of the 320 g EDELWEISS detector GGA1. The 
exposure of the fiducial volume corresponds to 8.6 kg-d. Also plotted as full 
lines are the ±1.645cr bands (90% efficiency) for photons and for nuclear 
recoils. The 99.9% efficiency region for photons is also shown (dotted line). 
The hyperbolic dashed curve corresponds to 3.5 keV ionization energy and 
the vertical dashed line to 20 keV recoil energy. 
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Figure 4: Spin-independent exclusion limits (dark solid curve) obtained by 
combining our 2000 data from Ref. || with the present data, for a total 
exposure of 11.7 kg-d. Dashed curve: previous EDELWEISS data || re- 
analyzed using the new definition of the upper bound of the recoil energy 
range (acceptance of 95%). Dash-dotted curve: present 2002 data. Closed 
contour: allowed region at 3a CL from the DAMAl-4 annual modulation 
data 0]. The full circle and triangle within this contour are defined in the 
text. 
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Figure 5: Combined EDELWEISS spin-independent exclusion limits (dark 
solid curve) compared with published limits from other experiments and the- 
oretical calculations. Dashed curve: Ge diode limit from IGEX ||. Dash- 
dotted curve: CDMS limit with statistical subtraction of the neutron back- 
ground @. Closed contour: allowed region at 3a CL from the DAMAl-4 
annual modulation data || . Two regions spanned by some of the supersym- 
metric model calculations of Refs. [H| are also shown. 
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